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The Convective Urea Process with Implicit Two-Dimensional 

Hydrodynamics 

J. Stein^’^ and J.C. Wheeler^’^ 

ABSTRACT 

Consideration of the role of the convective flnx in the thermodymics of the 
convective Urea nentrino loss process in degenerate, convective, qnasi-static, 
carbon-bnrning cores shows that the convective Urea process slows down the con¬ 
vective enrrent aronnd the Urea-shell, bnt, nnlike the “thermal” Urea process, 
does not reduce the entropy or temperature for a given convective volume. Here 
we demonstrate these effects with two-dimensional numerical hydrodynamical 
calculations. These two-dimensional implicit hydrodynamics calculations invoke 
an artificial speeding up of the nuclear and weak rates. They should thus be 
regarded as indicative, but still qualitative. We hnd that, compared to a case 
with no Urca-active nuclei, the case with Urea effects leads to a higher entropy 
in the convective core because the energy released by nuclear burning is conhned 
to a smaller volume by the effective boundary at the Urea shell. All else being 
equal, this will tend to accelerate the progression to dynamical runaway. We 
discuss the open issues regarding the impact of the convective Urea process on 
the evolution to the “smoldering phase” and then to dynamical runaway. 

Subject headings: Physical processes: convection - hydrodynamics - nuclear re¬ 
actions - Stars: interiors - supernovae 


1. Introduction 

The nature of degenerate carbon burning prior to runaway in Type la supernovae is a 
long-standing problem that has become more acute with the need to understand the variation 
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in the properties of the supernovae that might affect their application to the measurement 
of cosmological parameters. In this context, it is important to determine the effect of the 
convective Urea process because it may have direct observational implications for Type la 
supernovae. The convective Urea process can affect the density at which carbon undergoes 
dynamic runaway, subsequent electron capture and neutronizing reactions in the explosion, 
the size and gradient in the carbon/oxygen ratio at runaway, and the spectrum of convective 
motions and hence subsequent Rayleigh-Taylor unstable dynamic burning. 

The essence of the convective Urea process was hrst presented by Paezyhski (1972) 
whereby the convective circulation driven by carbon burning in degenerate white dwarfs will 
yield hrst electron capture and then /d-decay of susceptible nuclei. This yields no net change 
in composition, but a loss of neutrinos (or antineutrinos) at each step of the cycle along with 
their attendant energy. Paezyhski argued that this cycle would cool the star and postpone 
dynamical runaway. Previously, it had been pointed out by Bisnovatyi-Kogan & Seidov 
(1970) that non-equilibrium weak interactions as involved in this process would yield local 
heating. Bruenn (1973) independently pointed out that each weak interaction added heat 
to the system, despite the loss of the neutrino. As electrons are captured below the Fermi 
sea, another electron will drop from the Fermi surface to hll the “hole,” resulting in heat. 
In the other portion of the cycle, /5-decay will deposit electrons with excess thermal energy 
above the Fermi sea. Paezyhski (1974) acknowledged the arguments of Bisnovatyi-Kogan 
& Seidov and of Bruenn and pointed out that the energy flow between carbon burning and 
Urea neutrinos depends on the interaction between the convective and Urea processes. He 
noted that evolutionary calculations by Couch & Arnett (1973) and by Ergma & Paezyhski 
(1974) were done with the implicit assumption that Urea neutrinos directly remove heat 
from the stellar interior. In this paper Paezyhski also noted that the convective Urea process 
will lead to a gradient in the electron abundance. 

The physics of the convective Urea process was analyzed by Barkat & Wheeler (1990; 
hereafter BW), who argued that convective currents of composition play a critical role in 
the quasi-steady state thermodynamics of the convective Urea process, so that cooling terms 
associated with the convective flow exactly cancel the microscopic heating terms described 
by Bruenn and lead to net cooling. Their results seemed to be consistent with the results of 
the careful numerical work of Iben (1978a,b, 1982). 

The conclusions of BW were called into question by Mochkovitch (1996). Mochkovitch 
presented a general thermodynamic analysis and concluded that the convective Urea process 
must heat the star. These issues were addressed again by Stein, Barkat & Wheeler (1999: 
hereafter SBW) who concluded, in agreement with Mochkovitch, that while the neutrinos 
associated with the convective Urea process carry away energy, the entropy, and the temper- 
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ature, cannot decline. BW inclnded a “work” term that effectively removed energy from the 
total energy bndget that could only have come from the kinetic energy, which must remain 
positive. The loss rate by the Urea process cannot exceed the rate of generation of kinetic 
energy. Rather, the convective Urea neutrino losses slow the convective currents (SBW, 
Bisnovatyi-Kogan 2001; Lesaffre, Podsiadlowsi & Tout 2005). 

Here we present calculations demonstrating the role of the convective Urea process in 
two-dimensional simulations and illustrate that the conclusions of SBW were correct. In 
particular, we examine the role of convective currents, buoyancy and mixing in creating and 
limiting the kinetic energy. 

We present the outline of the input physics in §2 and the numerical technique in §3. 
The results are presented in §4 and conclusions are presented in §5. 


2. Description of Physics and Models 
2.1. Initial setting 

We constructed a Carbon-Oxygen core (1.38M©) with equal abundance of C and O by 
mass and with pc = 3.5 x 10® g cm“^ and Tc = 3 x 10® K in which the inner region {0.79Mq) 
had constant entropy, and the outer region was isothermal. For the ^®Ne/^®Na Urea pair, 
the Urea shell is at a density of 1.7 x 10® g cm“®. In these calculations, the Urea shell falls 
in the middle of the constant-entropy region, and divided that region into an inner zone 
(O. 3 IM 0 ) and an outer zone. This setup was necessary to distinguish the effects of the Urea 
process on the convection from the effect of the steep entropy gradient in the isothermal 
part. The inner zone contained a small concentration (0.00004 — 0.0004%) of “mother” Urea 
nuclei (^®Ne) and the outer zone with the rest of the star contained the same concentration 
of “daughter” Urea nuclei (^®Na ). This division between the inner and outer zones is stable 
against convection, even without the Urea process, as long as the entropy of the outer zone 
is equal to (or even very slightly less than) the entropy of the inner zone. In order to see 
the stabilizing effect of the Urea process on the convection we had to create a marginally 
stable conhguration. We achieved this by slightly reducing the entropy of the outer zone, 
and creating a slightly unstable configuration. This instability triggered a violent convection, 
which settled down after a short time, leaving the simulation in the desired configuration. 
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2.2. Reaction rates 

For the Carbon burning we used a simplified reaction rate which turns Carbon into 
Oxygen: 

^ = -0.5775X2 (p9)2-79(T9)22 

where pg is the density in units of 10 ® g cm“^ and T 9 is the temperature in units of 10 ® 
K. The q-value for the reaction is g = 2.857 x 10^^ erg gm“^. The Urea data for density 
thresholds and rates were taken from Ergma and Paezyhski (1974). In order to distinguish 
between the convective Urea process and the cooling thermal Urea process, we neglect the 
thermal Urea process. Other neutrino losses were taken from Itoh and Kohyama (1983). In 
our simulations we multiplied all the rates by factors 2000 or 20000 in order to speed up 
the process. Even so, the simulations required hundreds of thousands time-steps to produce 
~ 10 "^ seconds of simulation “star time.” 


3. Numerical Procedure 

The goal of this program is to simulate convection in a marginally stable carbon-oxygen 
core, where the sound-speed is high and the material velocity is low. Two-dimensional sim¬ 
ulation of subsonic convection in the central part of a white-dwarf presents some major 
stability problems: 

1. Having nearly the Chandrasekar mass, the star is only marginally stable to small pertur¬ 
bations. 

2. Convection occurs in a region with huge density and temperature gradients, but almost 
constant entropy. The convection is driven by the tiny entropy gradients. Small density or 
temperature fluctuations caused by numerical inaccuracies may create large entropy gradi¬ 
ents, destabilizing the convection. 

3. The velocities are subsonic, so V-(pu) is almost zero, in spite of the large density gradients. 

These issues consequently require a very stable numerical procedure. In order to ad¬ 
dress these stability issues, we therefore incorporated the following features into our program 
entitled DWARF: 

1. We use an implicit Eulerian scheme with second-order donor. 

2. In order to avoid the “chessboard black-white” decoupling of the cells in a usual stag¬ 
gered mesh, where density, energy and pressure are defined inside the cells and the velocity 
is defined on the corners, we define the velocities on the edge-centers of the 2D cells. 

3. Since the pressure differences between adjacent cells due to density differences may be 
larger by several orders of magnitude than the pressure difference due to the entropy differ- 



ences, we take special care that the advection terms will not create spurious large entropy 
gradients. 


3.1. Hydrodynamic equations 


The hydrodynamic equations are: 


Acceleration: 
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Internal energy: 
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where u is the velocity, p is the density, e is the energy per unit mass, G is the gravitational 
force and q is the energy source per unit mass per unit time. We use spherical coordinates, 
r, 6, 0, and assume that the gravity has spherical symmetry. All other quantities are treated 
in cylindrical symmetry. 

In these coordinates and under these assumptions our equations become: 


Acceleration: 
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We set rigid boundaries in angle at 6 = 45° and at 135° and in mass at an inner boundary 
of 10“^ Mq and an outer boundary of 1.38 Mq. The two angle and inner mass boundaries 
were chosen to avoid the singularity of spherical coordinates. The outer mass boundary was 
adopted for stability reasons. Despite the latter rigid outer boundary, the radial relaxation 
(expansion) of the convective region was not affected. 
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3.2. Features of the Program DWARF 

The DWARF code incorporates a rectangular grid in polar coordinates. It uses a 
staggered-mesh with implicit-pressure scheme, in which only the pressure is treated implicitly 
(See Livne 1993 and references therein), is adequate for subsonic flow in an otherwise stable 
medium. In our particular case, this model shows strong Chessboard black-white decoupling 
of the cells, causing diagonal instability. This decoupling is caused by the fact that in this 
model is nearly zero, when (/c, /) = (i, j ± 1) or (/c, /) = (i ± 1, j). To avoid this instabil¬ 
ity we chose a different staggered-mesh, in which the velocities and accelerations are defined 
at edge-centers. Other physical quantities are defined at the cell-centers. For the implicit- 
pressure scheme, (see Livne 1993 and references therein), only the 

matrix is inverted; however, this derivative accounts for both acceleration and advection. 

To compute the acceleration, we divide the acceleration into two parts: 

1- Direct accelerations, ~ fl'r — and are treated implicitly. 

2- Centrifugal (^u^) and Coriolis {—^UrUe) accelerations are treated explicitly, using the 
values at the end of the time-step. 

Stability requires that the advection be computed by a donor scheme; however, a donor 
scheme is difficult to treat implicitly, because the dependence of the advection on the velocity 
has discontinuous derivatives. Therefore we divided the advection into two parts: 

1- Using centered values (no donor) in the implicit part. 

2- Correction to a second order donor scheme is done explicitly. 

Choosing the “centered” donor can be done in several ways. Since the convection is driven 
by entropy gradients, we defined “centered” to be average density, entropy and composition 
of the two adjacent cells. The energy and temperature were calculated from these values by 
using the equation of state. 


4. Results 

We present the results of 5 simulations:^ 

S-4-2E4-0: (Xjve = 4 x 10“^), reaction-rates multiplied by 20000, Urea-rate = 0. 

S-4-2E4: (X^re = 4x 10“"^), reaction-rates multiplied by 20000, Urea-rate multiplied by 2000. 
S-1-2E4: (Xate = 1 X 10“"^), reaction-rates multiplied by 20000, Urea-rate multiplied by 2000. 
S-0.4-2E4: (X^re = 4 x 10“^), reaction-rates multiplied by 20000, Urea-rate multiplied by 
2000 . 


^computational models HN, FH, GG, EG and FI, respectively 
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S-4-2E3: (Xate = 4 x 10“^), reaction-rates multiplied by 2000, Urea-rate multiplied by 2000. 
When we multiplied the Urea-rate by another factor of 10 (to 20,000), the results did not 
change signihcantly, although the Urea-gradients became steeper, spreading over fewer cells. 
We thus left Urea rate multiplication factor at 2000 for all models where Urea effects were 
incorporated. 

All simulations were computed a long time after the convection settled down. Sim¬ 
ulations S-4-2E4, S-1-2E4 and S-0.4-2E4 (varying the Urea abundance) ran up to 6000 
simulation-time seconds. Simulation S-4-2E4-0 (with no Urea) ran up to 10,000 simulation¬ 
time seconds and simulation S-4-2E3 (with reduced nuclear rate acceleration), which was 
slower than the rest, ran up to 20,000 simulation-time seconds. The velocity scale in all the 
2D hgures (3-6 below) is normalized so that a single horizontal tic mark corresponds to a 
velocity of 10 km s“^. For the online color hgures the velocity is as follows: black arrows: 
below 10^ cm s“^; red arrows: 10^ - 2 x 10^ cm s“^; green arrows: 2 x 10^ - 4 x 10^ cm s“^; 
blue arrows: 4 x 10^ - 8 x 10^ cm s“h 

The energy input by nuclear reactions is converted in part to the kinetic energy of 
the convective circulation. The fraction of the input energy that ends up in kinetic energy 
depends on the dissipation to Urea “viscosity,” to turbulent viscosity and, for simulations, 
into numerical viscosity. For reference, the total energy input rate for simulation S-4-2E4 is 
2 X 10^^ erg s“^. The Urea neutrino loss rate is about 10^^ erg s“^. The total Urea dissipation 
rate including heating is about 1.3 x 10*^^ erg s“^, about 6.5 percent of the energy input. As 
shown in SBW, the rate of conversion of thermal energy to kinetic energy up to the Urea 
shell (T ~ 1.9 X 10® K) is less than (Tc — Tu)/Tc ~ 37 percent, where Tc is the central 
temperature and Tu is the temperature at the Urea shell. This implies that the Urea process 
as computed here makes a significant contribution to the kinetic energy dissipation in terms 
of the estimated maximum dissipation. 

Figure 1 shows the kinetic energy of all the simulations. In simulations S-4-2E4, S-1-2E4, 
and S-0.4-2E4 in which the nuclear rates were all multiplied by 20,000, but the abundance of 
Urea nuclei varied by a factor of 10, the kinetic energy associated with the convection settled 
rapidly into a quasi-steady state after recovering from the initial instability. In simulation 
S-4-2E4-0, which had no Urea effect, the convection spread to £11 the whole inner, constant 
entropy region, and so the approach to steady state took longer. Simulation S-4-2E3, with 
the nuclear rate enhanced by only a factor of 2000, took much longer to reach steady state. 
Simulations S-4-2E4, S-1-2E4, and S-0.4-2E4 converged to more or less the same kinetic 
energy, showing that in these simulations the kinetic energy does not depend sensitively on 
the abundance of the Urea nuclei, as long as the Urea physics is active. Simulation S-4-2E4-0, 
which lost no kinetic energy to the Urea effect, reached a somewhat higher value of kinetic 



energy. Simulation S-4-2E3, with a nuclear rate 10 times less asymptoted to a value around 
4 times smaller than simulations S-4-2E4, S-1-2E4, and S-0.4-2E4. This gives some notion 
of the dependence of the convective kinetic energy on the rate of nuclear energy input. 

Figure 2 shows the total Urea neutrino loss rate of all the simulations. Although there 
are large fluctuations, the general trend is that a higher abundance of Urea nuclei leads to 
sma/Zer neutrino losses. Specihcally, simulation S-4-2E4, which had the highest concentration 
of Urea atoms, lost less energy to Urea neutrinos by a factor of 2 to 3 than simulations S-1- 
2E4 and S-0.4-2E4. This is presumably related to the fact that the higher ^^Ne abundance 
restricted the extension of the convection beyond the Urea shell and the attendant mixing 
outward of carbon-burning products compared to simulations with small abundances of 
Urca-active nuclei (see Figures 3 below). On the other hand, we cannot rule out that this 
difference is affected by our numerical resolution. Simulation S-4-2E3, which had the smallest 
nuclear rate enhancement lost less energy to the Urea neutrinos than simulation S-4-2E4, by 
a factor of about 5, which is very similar to the ratio of a factor of 4 between the asymptotic 
kinetic energies of the two simulations. This suggests that the Urea neutrino loss rates are 
roughly proportional to the kinetic energy in the regime where Urea losses control the extent 
of convection (see the Appendix for a discussion of conditions where strong convection may 
overwhelm the Urea process). 

The question this work was designed to address is: “What does the Urea effect do?” 
The answer to this question, our main result, is presented in Figure 3. These hgures show 
the percentage of burned carbon (more accurately: 0.5 — Xc), for simulations S-4-2E4-0, 
S-4-2E4, S-1-2E4, and S-0.4-2E4 respectively. The carbon burns near the center. When no 
Urea effect is present - the convective currents spread the carbon-burning products evenly 
over the whole inner region (Figure 3A, simulation S-4-2E4-0). In Figure 3B, simulation S-4- 
2E4, only a small percentage of the carbon-burning products sneaks through the Urea-shell. 
This percentage increased as Xjve was decreased (Figures 3C and 3D). 

Figure 4 shows the entropy in simulations S-4-2E4-0 and S-4-2E4, respectively. Note 
the decrease of the entropy beyond the Urea-shell in simulation S-4-2E4, which compensates 
for the increase in the number of electrons per atom. 

Figure 5 shows the abundance of Urca-active Ne^^ atoms in the inner region of sim¬ 
ulations S-4-2E4-0 and S-4-2E4, respectively. As expected, in simulation S-4-2E4-0 the 
abundance is homogeneous, while in simulation S-4-2E4, the Ne^^ abundance is 0.0004 in¬ 
side the inner zone within the Urea shell, and zero in the outer zone beyond the Urea shell, 
except for a narrow transition zone near the Urea-shell. Figure 6 shows the Urea neutrino 
loss-rate in simulation S-4-2E4. The whole Urea activity occurs near the Urea-shell, with 
the exception of some boundary-related activity. In all the 2D hgures, the overshootings at 



the left and right boundaries are numerical artifacts. 


5. Conclusions 

These two-dimensional calculations make no assumptions about the convective Urea 
process per se; they employ the relevant weak interaction microphysics and then follow the 
affects of the convective flow with no extra assumptions about thermal or chemical equi¬ 
librium. Work done by convection and associated neutrino losses are computed straightfor¬ 
wardly. Because the computations have been done with artihcially accelerated thermonuclear 
and weak interaction rates, they cannot be taken literally; some processes like turbulent dis¬ 
sipation will not scale linearly with these rates. It would, of course, be desirable to do these 
simulations in three, rather than two dimensions. We expect the qualitative results to be 
the same, but this must eventually be checked. Our results should therefore be judged for 
their qualitative value. 

With these caveats, we conclude here: a) In agreement with Mochkovitch (1996) and 
SBW and in disagreement with the earlier work in BW, that the convective Urea process 
cannot result in a net decrease in entropy, and hence in temperature, for a constant or 
increasing density; b) that the Urea process controls the convective flow and affects the 
kinetic energy even as the convective motions drive the Urea process; c) that the convective 
Urea process will, at least temporarily, limit the expansion of the convective zone beyond 
the Urea shell. Bisnovatyi-Kogan (2001) also noted that the convective Urea process must 
slow the speed of convection. Despite their qualitative nature, we feel these simulations have 
established that the convective Urea process is fed by the positive definite kinetic energy and 
cannot cool the star, but can and will affect the nature of the convective flow. 

Over the last couple of decades most work on modeling of the progenitors of Type la 
supernovae has simply ignored the effects of the convective Urea process. It remains to 
be seen whether the convective Urea process will have a signiheant effect on the ultimate 
dynamical runaway. The convective Urea process by itself will not lead to dynamical runaway 
at higher density because it cools the star, as hypothesized in some earlier work. The 
convective Urea process may shorten the “smoldering” phase (Hoflich & Stein 2002) between 
carbon ignition, dehned as the epoch when carbon burning exceeds neutrino losses, and 
dynamic runaway. Compared to a model without the effect of the convective Urea process, 
the convective core, and hence the heat liberated by the carbon burning, is conhned to a 
smaller volume. All else being equal, a model that invokes the convective Urea process should 
proceed to dynamical runaway more quickly than one without the convective Urea process. 
The net result might be to indirectly lead to dynamic runaway at lower density with effects 
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on the explosion dynamics, nucleosynthesis, and light curve. 

These issues will depend sensitively on whether there is an active Urea shell and that, 
in turn, depends on the density of carbon ignition. While that issue was once actively 
explored, it has received little commentary in recent literature on the topic of degenerate 
carbon burning. This is in part due to advances in understanding of weak interaction rates 
(see Brachwitz et al. 2000, and references therein) and partly due to settling on a standard 
set of assumptions for the treatment of effects such as plasmon neutrino losses and strong 
screening effects (that may or may not be correct). 

Calculations of non-rotating accreting white dwarf evolution tend to give carbon ignition 
at densities exceeding 4x 10 ® gm cm“^ for modest accretion rates, ~ 10 “® M 0 yr“^, depending 
on the mass of the original white dwarf and the temperature at the beginning of the accretion 
phase (Bravo et al. 1996; Hoflich, Nomoto, Umeda & Wheeler 2000). Lesaffre et al. (2006) 
estimate the distribution of central densities for single-degenerate progenitors and hnd a 
broad distribution with a minimum ignition density of ~ 2 x 10 ® g cm“^, above the threshold 
density for the Ne/Na Urea pair. Such conditions will inevitably be accompanied by one or 
more convective Urea shells. These conditions might be expected to lead to nova outbursts 
that prevent growth to runaway. If growth of the white dwarf continues, the result could 
be accretion-induced collapse (Saio & Nomoto 1998). It will be interesting to investigate 
whether or not conhnement of convection to Urea shells and the subsequent tendency to heat 
a smaller volume more rapidly would affect this conclusion. Only for rather high accretion 
rates, approaching 10 “® M 0 yr“^, will the central density at carbon ignition be less than the 
threshold density for the Ne/Na convective Urea cycle so that the convective Urea process 
might not occur at all. For these high rates, there is some danger of excessively polluting 
the circumstellar environment with a hydrogen-rich wind. 

Observational constraints on Type la supernovae tend to put the central density at igni¬ 
tion in a range where the Urea process might occur, independent of how the star (or model) 
got there. Plausible calculations and their dynamical outcomes that seem to reproduce ob¬ 
servations of Type la supernovae (Hoflich, Nomoto, Umeda & Wheeler 2000; Dominguez, 
Hoflich & Straniero 2001) give ignition densities of about 1.5 —2.5 x 10® gm cm“^, right where 
the Ne/Na convective Urea “on” switch might lie. This means that whether a model (and 
hence a real star, in principle) does or does not trigger the convective Urea process with its 
potential effect on convective velocity fields, extent of convection, amount of carbon burning 
between ignition and runaway, and the C/0 gradient at runaway, may depend sensitively 
on such recently neglected issues as electron screening and plasmon neutrino loss rates. We 
also note that the white dwarfs that explode as Type la supernovae will almost surely be 
rotating. The rotational state will also affect the carbon ignition density (Piersanti et al. 



2003a,b) and hence the presence of convective Urea effects. 

Another important issue is the abundance of a given Urca-active pair. We have shown 
that the kinetic energy of convective motion is rather independent of the abundance of the 
Urea nuclei, but that the mixing of carbon and its burning products does depend on the Urea 
abudance. We also hnd that neutrino losses are reduced (within the range of parameters 
explored here) for higher abundance of Urea nuclei, although we cannot rule out that this 
is an effect of our numerical resolution. The abundance of Urca-active nuclei are a function 
explosive nucleosysthesis in prior generations of stars and so the role of the convective Urea 
process might be a function of redshift. Another possible effect on abundances is evolution 
in the core. If burning prior to runaway should produce more Urea nuclei, then the Urea 
process could be enhanced. 

Other interesting issues are related to the location of the Urea shell (or shells) and the 
speed of convection. How does the Urea process work if the central density is only slightly 
above the Urea threshold so that the Urea shell is quite small, smaller, perhaps, than a 
convective scale length? In general, how fast is the convection? Sufficiently rapid convection 
will sweep mothers back down to regions where they are stable before they can decay. We 
note that when carbon ignition first occurs, the convection is likely to start slowly, giving the 
Urea process ample time to work, as long as the ignition density exceeds the Urea threshold. 

The key question is, even if there is an active convective Urea stage in the evolution 
toward dynamic runaway, does this make any substantial difference in the nature of that 
runaway and the subsequent dynamics? If the Urea process limits the radial extent of the 
convection, somewhat more carbon might be consumed within the Urea shell than beyond. 
On the other hand, the total amount of carbon consumed is likely to be small and if the 
limitation on the extent of convection breaks down before dynamic runaway, then there will 
be a homognenization process as the convection sweeps outward and dredges carbon inward. 
It seems unlikely that an effect of the convective Urea process on the abundance, per se, will 
be signiheant. 

A potentially more interesting effect is on the velocity structure of the convective mo¬ 
tion. Hoflich & Stein (2002) investigated the hnal 3 hours before dynamic runaway with 
implicit two-dimensional calculations. They showed that the “smoldering” phase of nuclear 
burning that immediately preceeds dynamical runaway can produce convective motions with 
velocities of order 40 - 120 km s“^. There were two implications of this work. One was that 
dynamical runaway was triggered at a single region near to, but definately away from, the 
center of the star by compression driven by converging convective plumes. The second im¬ 
plication was that the convective speeds were comparable to, or even exceeded, the velocity 
of the burning front in the early, slow, phase of Rayleigh-Taylor unstable turbulent deflagra- 
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tion. This means that the convective velocity held can substantially “prepare” the conditions 
prior to dynamical runaway, can trigger the dynamical runaway, and that in the early phases 
of dynamical runaway it is still the convective velocities, not the burning front velocity, that 
determine the speed of propagation of the burning front. Recent calculations by Kuhlen 
et ah (2005) have addressed aspects of the late smoldering phase in three-dimensional cal¬ 
culations. They have not followed the burning completely to dynamic runaway, but hnd 
qualitatively similar results over the limited time frame they investigated: convective hows 
of order 50 - 130 km s“^ and hot spots near the center where dynamic runaway might occur. 
The outcome of dynamical ignition may depend sensitively on the number and location of 
the runaway spots leading in some cases to single plumes of burning and in others to more 
homogeneous outward propagation of the early dehagration burning (Livne, Asida & Hohich 
2005). 

The convective Urea process could thus ahect the nature of the convective velocity held 
and hence the smoldering phase, the location of dynamic runaway, and the early stages of 
dynamical runaway. The current calculations allow us to raise such issues, but because of the 
artiheial speed up of the nuclear and weak interaction rates, we cannot address them directly 
here. The speed up of the nuclear rates by a factor of 20,000 will result in a speedup of the 
convective velocities. In the current models, the convective motions are of order 10 km s“^. 
The question is what happens in a realistic model. Will the presence of the convective Urea 
shell continue to conhne the convection and will the convective velocities in the smoldering 
phase that precedes dynamical runaway be restricted and hence less able to affect the early 
stages of deflagration? Or will the convection speed up and somehow “break out” of the Urea 
shell yielding results similar to those of Hoflich & Stein (2002) who neglected the convective 
Urea process. 

In this context, the important issue is whether the convection becomes strong enough 
within an active Urea zone to overwhelm the Urea control and thus spread the convection 
throughout its “natural” regime as if there were no Urea processes. The convective Urea 
process might control the extent of the convection through much of the smoldering phase 
only to be overwhelmed in the hnal few hours or minutes before dynamic runaway; sufficient 
time to establish a convective structure that lost all memory of the convective Urea phase. 

That this breakdown of the convective Urea process could occur can be seen qualitatively 
by considering that for sufficiently strong convection, convective flows may push substantially 
beyond the Urea shell to the point where all the mothers have decayed to daughters and 
there will be no more neutrino losses regardless of how far the convection thus extends. We 
have shown in Figures 1 and 2 that the convective kinetic energy and neutrino loss rates 
increases with the nuclear reaction rates, but, at least for these simulations, at a rate that is 
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slower than linear. These simulations also suggest that within the range we have explored, 
the neutrino loss rates scale roughly linearly with the convective energy. These conclusions 
hold as long as the Urea shell can control the convection. We sketch in the Appendix 
an argument that suggests that sufficiently strong convection can cause the convection to 
“break out” and overwhelm the limits of the convective Urea process. In principle, our use 
of accelerated nuclear reaction rates might correspond to a later phase in the smoldering 
process, closer to dynamical runaway, although it is not clear this is true. Based on our 
simulations, we estimate in the Appendix that convective velocities substantially in excess 
of 20 to 70 km s“^, depending on the abundance of the Urca-active nuclei, are required for 
“breakout.” Interestingly, this is similar to the velocity range found in the studies of Hoflich 
& Stein (2002) and of Kuhlen et al. (2005) for the late smoldering phase. This suggests 
that breakout might, or might not, occur prior to dynamical runaway. We conclude that 
this issue is open and in need of thorough future study. While there is no obvious reason to 
think the convective Urea process will affect subsequent dynamic phases, we currently have 
no grounds to rule out that it may do so. 

Although the convective Urea process has been long neglected, the current calculations 
suggest that it needs to be in the suite of actively-considered physical phenomena as we 
attempt to achieve a rehned level of understanding of Type la supernovae. We have iden- 
tihed a number of important issues regarding the convective Urea process that need to be 
explored in multidimensional calculations that investigate the evolution from carbon ignition 
to dynamical runaway. A key issue is the behavior of the convective Urea process without 
the artiheal enhancement in rates we have employed and in three dimensions. 

We are grateful to Peter Hoflich for discussion of related topics of degenerate convective 
carbon burning. This research was supported in part by NSF Grant AST-0098644, 


A. An Estimate of Urea Stopping Power 

In order to get a rough estimate of the stopping power of the Urea effect, let us assume 
that one blob of unit mass of material containing concentration of mothers, crosses the 
Urea-shell outward, its velocity u is perpendicular to the Urea-shell. Let us further assume 
that the Urea effect is local, i.e., when mothers become daughters in the blob, the blob 
loses kinetic energy equivalent to the neutrino -|- heat energy released in the process. This 
estimate would be straight forward except that the energy lost depends on the location in 
which it was released: the further from the Urea-shell, the more energy is lost to neutrinos 
plus /5-decay heating. 
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Let X be the distance of the blob from the Urea-shell, then, nnder the above assnmptions, 
the motion of the blob is governed by the following eqnations: 


dx 

dt 


(Al) 


^ (A2) 

dt 

where C^vi^x^ is the rate of captnre per atom with Cy being a nnmerical coefficient (nnits 
erg-3 g-i^ characteristic of a given Urea pair, rj = |^^|, at a; = 0, Na is Avogadro’s nnmber 
and A is the atomic nnmber of the Urea atoms. 


If the blob is not stopped by the Urea process, bnt loses a fraction of its kinetic energy, 
then Xm goes to 0 for large x where all /5-decay has gone to completion. Then has a 

maximnm at some point Xc- Let us assume that the whole reaction occurs near that point. 
At the maximum we have: Xc = The energy released at the maximum is: 


TA _Na^ ( 1 _Na^ ,^ur] 1 

-^C — ^ XmV^c — ^ XmV ^ ~ Q )' 


(A3) 


Therefore, in order to continue moving outward, the blob should have kinetic energy signif¬ 
icantly larger than E^, or 

/V A 1 A^iin 1 

(A4) 


.lASur] 1 
n > (— Xj 2(_^)8 


or 


u > 


ASVaNa . 4 . 


(A5) 


Note that this criterion depends only weakly on the Urea rate and is only mildly sensitive 
to the Urea abudance. 


In our case z/ = 2000 x 3.3 x 10^^ = 6.6 x 10^® erg“^ s“^ (Ergma & Paezynski 1974), 
f] = 0.86 X 10“^^, and A = 23, so when X^ = 4 x 10“^ then u must be significantly greater 
than 70 km s“^. When X^ = 4 x 10“^, u must be signihcantly greater than 19 km s“^. 
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Fig. 1.— Plot of convective kinetic energy versus simulation time, for simulations S-4-2E4, 
S-1-2E4, S-0.4-2E4, S-4-2E3, and S-4-2E4-0. Simulations S-4-2E4, S-1-2E4, S-0.4-2E4, and 
S-4-2E4-0 had the nuclear rates multiplied by 2 x 10^ and simulation S-4-2E3, had the nuclear 
rate multiplied by 2 x 10^. 
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20000 


Fig. 2.— Plot of Urea neutrino loss versus simulation time, of simulations S-4-2E4, S-1- 
2E4, S-0.4-2E4 and S-4-2E3. Simulations S-4-2E4, S-1-2E4, S-0.4-2E4 had the nuclear rate 
multiplied by 2 x 10"^ and simulation S-4-2E3 had the nuclear rate multiplied by 2 x 10^. 
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Fig. 3.— A) (upper left) Plot of the distribution of carbon-burning products, velocity vectors, 
the Urea shell (heavy black line around R = 400 Km) and the boundary separating the 
inner region of constant entropy from the outer isothermal region (thinner black line around 
670 Km) for simulation S-4-2E4-0 with no Urea process. The carbon-burning products are 
distributed more or less evenly in the whole inner region. B) (upper right) Same as A but 
for simulation S-4-2E4 with the Urea process turned on . The carbon-burning products are 
almost entirely conhned to the inner zone (within the Urea-shell). Strong convection is seen 
in the inner zone, and a separate convection region exists in the outer zone between the 
Urea shell and the boundary of the isothermal region, probably triggered by the overshoot 
from the inner zone. C) (lower left) same as B) but for simulation S-1-2E4 with smaller 
abundance of Urca-active nuclei. The carbon-burning products are less conhned to the inner 
zone (within the Urea-shell) than in B. The velocities near the Urea-shell are almost parallel 
to the shell. D) (lower right) Same as B and C but for simulation S-0.4-2E4 with even 
smaller abundance of Urca-active nuclei. The carbon-burning products are less conhned to 
the inner zone (within the Urea-shell) than in C. Velocity vectors are normalized such that 
one horizontal tic mark represents 10 km s“^ (journal black and white version) or black 
arrows: less than 1 km s“^; red arrows: 1-2 km s“^; green arrows: 2-4 km s“^; blue 
arrows: 4-8 km s“^ (on-line color version). 
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Fig. 4.— A) (left panel) Plot of entropy, velocity vectors, Urea shell as in Figure 3 at 10,000 
seconds for simulation S-4-2E4-0 with no Urea effect. B) (right panel) same as A, but with 
Urea effect and time = 6000 sec for simulation S-4-2E4. 
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Fig. 5.— A) (left panel) Plot of /Na?'^ abundance, velocity vectors, Urea shell as in 
Figure 3 at 10,000 seconds for simulation S-4-2E4-0 with no Urea effect. B) (right panel) 
same as A, but with Urea effect and time = 6000 sec for simulation S-4-2E4. 
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Fig. 6.— Plot of Urea neutrino loss for simulation S-4-2E4. The whole Urea aetivity oeeurs 
near the Urea-shell. Areas adjaeent to the left and right boundaries are subjeet to numerieal 
artiheats and should be ignored. Veloeity veetors are represented as in Figure 3. 






